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^ ABSTRACT 
O 

<N NGC 3783 was observed for approximately 210 ks by Suzaku and in this time 

O showed significant spectral and flux variability at both short (20 ks) and long (100 ks) 

O time scales. The full observation is found to consist of approximately six "spectral 

OO periods" where the behaviour of the soft (0.3-1.0 keV) and hard (2-10 keV) bands are 

somewhat distinct. Using a variety of methods we find that the strong warm absorber 
present in this source does not change on these time scales, confirming that the broad- 
band variability is intrinsic to the central source. The time resolved difference- spectra 
are well modelled with an absorbed powerlaw below lOkeV, but show an additional 
hard excess at ^ 20 keV in the latter stages of the observation. This suggests that, 
O in addition to the variable powerlaw, there is a further variable component that varies 

^ with time but not monotonically with flux. We show that a likely interpretation is 

■ that this further component is associated with variations in the reflection fraction or 

^ possibly ionization state of the accretion disk a few gravitational radii from the black 

> hole. 

(N 
(N 

^ 1. Introduction 

^ Active Galactic Nuclei (AGN) are commonly found to display persistent X-ray spectral and 

flux variability in a variety of time-scales (e.g. |McHardy||1989[ [Green et al.|| 1993D. Studies of this 
broadband variability usually suggest that at least part is due to strong gravitational effects close to 
^ the central black hole. Time-resolved spectral analyses of such sources are thus an important tool 
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in the study of the inner regions of the accretion flow around supermassive black holes, providing 
information on both the geometry and physical evolution of the inner accretion disk and on intrinsic 
physical parameters such as black hole mass and spin. 

It has been shown for a number of Seyfert galaxies, most notably MCG-6-30-15 (see e.g. 



Fabian & Vaughan]|2003| ), that the spectrum above a few keV becomes harder at lower X-ray 
fluxes. At all flux levels, however, the X-ray spectrum of such systems usually consists of a pri- 
mary continuum which is accurately approximated as a powerlaw, together with a broad Fe-Ka 
line and cold reflection features ( |Nandra et al.||1997| ). In some cases it is also accompanied by 
highly ionized warm absorbers and a "soft excess" ( |Halpern|[l984[ |Reynolds]|1997[ ). The most 
common interpretation for the spectra of such AGN is that the primary X-ray emission irradiates 
the underlying optically thick material in the innermost part of the accretion disk, resulting in 
"reflection signatures" consisting of fluorescent and recombination emission lines as well as ab- 



sorption features. The most prominent signature is the broad, skewed Fe-Ka line (see e.g. Miller 



2007] for a recent review of relativistically broadened lines). Furthermore, the soft excess is a nat- 
ural consequence of reflection from the innermost regions around a black hole for a wide range of 
ionisation states. 



Recent spectral variability studies of Seyfert galaxies such as MCG-6-30-15 ( [Miniutti et al 
2007| ), 1H0707-495 ( [Fabian et al.||2009t |Zoghbi et al.|[20T0| ) and NGC 4051 ( jUttley et al.||2004 



Terashima et al.|[2008l ) have found different explanations for the origin of the variability - some- 



times in the same object. Amongst the possibilities are (i) the spectrum consists of a soft com- 
ponent (powerlaw) with constant spectral slope and variable flux, together with a constant hard- 
component associated with strong reflection from the innermost accretion disk (e.g MCG-6-30-15; 
Vaughan et al.||2003| [Taylor et al.||2003[ [Miniutti et al.||2007| ); (ii) the soft component consists of a 



pivoting powerlaw with the photon index increasing with the logarithm of the flux, together with a 
constant reflection component (e.g. NGC 4051; [Uttley et al.[[2004[ ); and (iii) a pivoting powerlaw 



together with a constant, partially-covered, neutral (cold) reflection (e.g. NGC 4051; [Terashima 
et al.||2008| ). Suzaku provides an ideal opportunity to discern between these various distinct inter 



pretations due to its broad-band spectral coverage and good signal-to-noise in the full 0.5-50 keV 
range as well as large effective area below ^10 keV. 

NGC 3783 is a bright, nearby Seyfert 1 galaxy exhibiting prominent broad emission lines 
and strong X-ray absorption features. Based on a total of approximately 900 ks of Chandra ob- 
servation, various authors concluded that this source contains three zones of ionized absorption 
covering a large range of ionization ( Kaspi et ar[|2002[ jNetzer et"ar[[2003[ jKrongold et al.[|2003t 
[20051). Analyses of the grating data showed little, if any, evidence for changes in the optical depths 



of the absorbers. Similar conclusions have also been made using XMM-Newton RGS data ( [Blustin 
et al.[[2002| ). We observed NGC 3783 for approximately 340 ks with Suzaku as part of the Suzaku 
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AGN Spin Survey Key Project. A detailed analysis of the time averaged spectrum is presented in 
Brenneman et al.| ( [2011| ) (hereafter Paper 1). For the first time in NGC 3783, the dimensionless 



spin parameter for the central black hole was estimated to be > 0.88 at the 99 per cent level of 
confidence. This was based on the self-consistent modelling of the full reflection features, whilst 
simultaneously accounting for absorption due to the warm absorber in the line of sight. In this pa- 
per we present an examination of the X-ray spectrum and spectral variability of NGC 3783 using 
this same long Suzaku observation with the goal of better understanding the accretion disk flow ge- 
ometry and evolution of this system. By usage of a variety of techniques we show that the spectral 
and flux variation seen in NGC 3783 are likely due to a combination of a varying powerlaw (with 
constant index) as well as a hard-reflection component. Furthermore, we find that the ionization 
and reflection fraction in the innermost region of the accretion disk are also likely to be varying. 

We start in §|2]with specifics of the observation and data reduction procedure. This is then 
followed by a detailed analysis of the spectral variability (§[3]) using a combination of flux-flux 
relations (§ |3.2| ), flux-resolved (§ |3.3| ) and time-resolved difference- spectra (§ |3.4| ). We apply the 



results of the previous sections to each individual time period in (§ |3.5.1| ). Our results are then 
summarised in §|4} 



2. Data reduction 

We observed NGC 3783 for a total of ^340 ks (yielding 210 ks of good "on source" expo- 
sure) with Suzaku ( [Mitsuda et al.||2007| ) as part of the Suzaku AGN Spin Survey Key Project. The 



observations were performed during the period 10-15 July 2009. The three operating detectors 
constituting the X-ray Imaging Spectrometer (XIS; |Koyama et al.|2007| ) were operated in the "nor- 



mal" clock and the 3x3 and 5x5 editing modes. The observation resulted in a total (co-added) 
exposure of approximately 420 and 210 ks for the front-illuminated (FI) and back-illuminated (BI) 
instruments, respectively. 

Using the latest HEASOFTv6.10 software package we processed the unfiltered event files 
for each of the XIS CCDs and editing modes operational, following the Suzaku Data Reduc- 
tion Guided We started by re-running the Suzaku pipeline with the latest calibration (2011-02- 
10), as well as the associated screening criteria files in order to create new cleaned event files. 
XSELECT was used to extract spectral products from these event files. The good time intervals pro- 
vided by the XIS team were employed in all cases to exclude telemetry saturations. Source events 
were extracted from a circular region of 200 pixels (^ 208^0 radius centred on the point source. 



http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/ 
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and background spectra from another region of the same size, devoid of any obvious contaminating 
emission. The script XISRESP^, which calls the tools XISRMFGEN and XISSIMARFGEN, was used 
with the "medium" input to obtain individual ancillary response files (arfs) and redistribution ma- 
trix files (rmfs). Finally, we combined the spectra and response files from the two front- illuminated 
instruments (XISO and XIS3) using the FTOOL ADDASCASPEC to increase signal-to-noise. The 
spectra, background and responses were then rebinned by a factor of eight to 512 channels using 
RBNPHA and RBNRMF . The FTOOL GRPPHA was used to give at least 100 counts per spectral 
bin. Throughout this paper BI and FI spectra are fitted in the 0.5-8.0 keV and 0.7-10 keV energy 
range respectively. Unless stated otherwise, the 1.7-2.3 keV energy is ignored due to the possible 
presence of uncalibrated instrumental features. All energy axes are in the observed frame. 

For the hard X-ray detector (HXD, [Tikahashi et al.|2007| ) we again reprocessed the unfiltered 
event files following the data reduction guide (only the PIN data are used in this analyses). Since 
the HXD is a coUimating rather than an imaging instrument, estimating the background requires 
individual consideration of the non X-ray instrumental background (NXB) and cosmic X-ray back- 
ground (CXB). The appropriate response file was downloaded^ together with the tuned (Model D) 
background. A common good-time interval was obtained with MGTIME which combines the good 
time of the event and background file, and XSELECT was used to extract spectral products. Dead 
time corrections were applied with HXDDTCOR, and the exposures of the NXB spectra were in- 
creased by a factor of ten, as instructed by the data reduction guide. The contribution from the 
CXB was simulated using the form of |Boldt| ( |1987| ), with the appropriate normalisation for the 
HXD nominal pointing. The NXB and CXB spectra were then combined using MATHPHA to give 
a total background spectrum, to which a 2 per cent systematic uncertainty was added. The PIN 
spectrum was grouped to have a minimum of 500 counts per energy bin to improve statistics, and 
again allow the use of minimization during spectral fitting. The PIN data reduction yielded a 
total source rate of (3.27 ±0.02) x 10"^ count s"^ which is 38.5 per cent of the total observed flux, 
with a good exposure time of 234 ks. We note here that since the main purpose of our analysis 
is to investigate the spectral variability of NGC 3783, all PIN spectra presented in this paper are 
constrained to be strictly simultaneous with their XIS counterpart. Furthermore, we restrict all our 
PIN analyses to the energy range 15.0^0.0 keV and fit it simultaneously with the XIS data by 
adding a normalization factor of 1.18"^ with respect to that of the FI spectrum, as recommended in 
the Suzaku Data Reduction Guide. 



^http://suzaku. gsfc.nasa.gov/docs/suzaku/analysis/xisresp 
^ http ://www. astro .isas . ac .jp/suzaku/analy sis/hxd/ 

"^See Paper 1 for a discussion of the effects in letting this cross normalization vary. 
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Fig. 1 . — Time-averaged data/model ratio to an absorbed powerlaw. XIS back and front illuminated 
spectra are shown in black and red respectively. The PIN data are shown in green. A powerlaw was 
fitted in the 3.0-4.0 and 7.5-40.0 keVenergy range, The residuals shows the various complexities 
including the presence of a strong warm absorber seen at soft energies together with a superposition 
of emission lines at ^ 6.4 keV. 

3. Data Analyses and Results 

Fig. [T] shows the time-averaged XIS/PIN spectra for the entire data set compared to a simple 
absorbed powerlaw model (A^h frozen at the value of 9.91 x 10^° cm"^ as found in Paper 1) fitted 
between 3.0-4.0 keV and 7.5-40.0 keV. The presence of a strong warm absorber is clearly seen 
below 2keV as is the presence of various emission features peaking at 6.4 keV. Detailed fits to 
the full Suzaku spectrum were presented in Paper 1. It was found that the time averaged spectrum 
was fitted well by a model consisting of a power-law continuum, blackbody-like soft excess, and 
X-ray reflection from both cold/neutral distant material and an ionized accretion disk; the observed 
emission from these components are also strongly affected by a three zone warm absorber. In order 
to adequately describe the time-averaged spectrum, a fraction (^15%) of the continuum emission 
scatters around or leaks through the warm absorbing gas. See Paper 1 (Table 1) for a detailed 
discussion of the time-average spectral model. 

In the following sections we will concentrate on the spectral variation of NGC 3783, initially 
using a number of model-independent techniques. 
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Fig. 2.— Top: 0.2-12.0 keV XIS light-curve in intervals of 1 ks. Middle: 15-40 keV PIN light- 
curve in intervals of 5 ks. Bottom: Normalised soft (0.3-1.0 keV, red), hard (2.0-10.0 keV, green) 
and soft/hard (colour, blue) count-rate as a function of time. In the lower panel these were nor- 
malised by dividing the difference between the count rate in each 20 ks time bin and mean count 



rate by the standard deviation. For § |3.4| the light-curve has been divided into six intervals (PI to 
P6 indicated by the dotted vertical lines. We also indicate intervals where both the hard, soft, and 
softness ratio are at distinct high and low levels. 



3.1. Colour changes 

It can be seen from the background- subtracted light curve in Fig. |2] (Top) that the X-ray 
continuum emission in NGC 3783 shows clear flux variability on a variety of time scales. The light 
curve shown above for the XIS data has bin sizes of 1 ks and is the result of the combined data 
for the three XIS instruments. The PIN light curve in the 15-40 keV range is shown in the middle 
panel with bin sizes of 5 ks. In the bottom panel we show the normalised 0.3-1.0 keV (red; soft) 
and 2.0-10.0 keV (green; hard) count rate together with the colour (blue), defined here as the ratio 
between soft and hard count rate. In order to better show the evolution of the colour, soft and hard 
energy band, we have spht the data into 20 ks time intervals. During the full observing campaign, 
NGC 3783 initially dropped in flux by a factor of ^ 1.8 in the first 100 ks and subsequently went 
back to an intermediate flux level during the latter ~ 200 ks. When comparing the evolution of the 
soft, hard and softness ratio, we see that there are approximately six distinct periods (P1-P6) which 
have been highlighted in Fig. [2} During the first and last two periods the softness ratio appears to 
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foUow both hard and soft normahsed rate variation. However, this is not the case in periods 3 
and 4 where the softness ratio is significantly below and above the two normalised flux bands, 
respectively. In order to characterise these spectral variations in a model-independent manner, we 
begin by applying the "flux-flux" analysis outlined in [Taylor et al.|(|2003|). 



3.2. Flux-flux analysis 



[Taylor etltTI ( |2003[ ) showed that for Seyfert galaxies, the flux measured in a specific energy 
band, F{E), can be expressed as the sum of constant {C^ and C\,) and variable (Fg and Fh) com- 
ponents, so that the flux in a hard energy band (subscript h), can be related to that of a soft band 
(subscript s), by the expression: 



(1) 



If the varying component has a constant spectral shape, the flux-flux relation will be linear 
(a = 1). If, however, the spectral variability is intrinsic to a single pivoting powerlaw with no 
contribution from any other component then plotting Fh versus will result in a simple powerlaw 
of the form F\, = kF^ where a can be used to estimate the pivot energy ( [Zdziarski et al.|2003| and 
Fig. 1 of [Taylor et al.|2003D . 

Figure [3] (Top) shows the hard (2.0-10.0 keV) versus soft (0.3-1.0 keV) count rates (flux- 
flux plot) for NGC 3783 from the 1 ks binned XIS light-curves. The scatter is partly due to poor 
statistics and is mostly removed when we bin the data. The red line is the best linear model which 
is summarised in Table [TJ Although visually the linear relationship appears correct, the fit is not 
statistically acceptable, with = 40.4/12. When we add a systematic error of 2 per cent to 



each XIS count (as per Noda et al. 2011), the resulting fit gives x^/^ = 11-5/12. This extra 
systematic is likely due to uncertainties in CCD contamination ( [Noda et aLl[2011| ). Replacing the 
linear relation with equation [T] results in only a mild improvement (x^/^ = 8.2/10; blue curve). 
This best fit cannot rule out weak pivoting nor does it statistically require it {a = 0.63 ± 0.16 with 
constant offsets on the hard and soft axes of C\, = 3.2 ± 0.2 and Cs = 0.54 ± 0.02 respectively; 



la). According to Figure 1 of [Taylor et al.[ ( [2003] ), a ranging from approximately 0.55 to 0.8 
implies a powerlaw with a pivot energy ranging from 30 to 300 keV. In Section [34] we show using 
time-resolved spectra, evidence for a distinct change in the spectrum between the first and second 
half of the observation. In order to investigate whether any change is apparent in the flux-flux 
plot between these two time intervals, we created flux-flux plots for the first 180ks (Pl-i-P2-hP3 
in Fig. [2]) and the latter 160 ks (P4-hP5+P6). These are shown Figure [3] (Bottom left and right 
respectively) and the fits are summarised in Table [TJ 



The total flux range probed in this observation is quite modest 1.8 between highest and 
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Fig. 3.— Flux-flux plots for NGC 3783. Hard (2.0-10.0 keV) versus soft (0.3-1.0 keV) count rate. 
Top panels: The points are from the light curve having 1 ks time bins. Bottom panels: Same as 
above but after averaging over the points and adding 2 per cent systematic to each XIS count. In 
all cases, the red line shows the best linear fit to the data and the blue curve represents the best fit 
described by Equation 1 . 

Table 1 : Summary of linear and powerlaw fits to the flux-flux plots shown in Fig [3j All errors are 
Icr. " The powerlaw is of the type shown in Equation [TJ 



Period 


Linear model 




Powerlaw model*^ 






Mill C 


X^/d.o.f 


h 

'"'pow 


a 


X^/d.o.f. 


Full range 


2.7 ±0.2 2.1 ±0.2 


11.5/12 


2.5 ±0.2 


0.63 ±0.16 


8.2/10 


0-180ks 


2.96 ±0.86 2.1 ±0.7 


9.9/6 


4.3 ±0.3 


0.55 ±0.22 


9.5/4 


180-340ks 


2.4 ±0.2 2.3 ±0.1 


7.3/13 


2.2 ±0.2 


1.05 ±0.07 


7.2/11 
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Fig. 4. — Binned band flux F(E2) versus 1.0-1.5 keV flux. The various lines shows the best 
linear-fit to a sample of energy range. By obtaining the constant offset for the various energies 
with respect to a standard range we can deduce the shape of the constant component (§ |3.2[ see 
Fig.|5]). 



lowest flux bins) compared with other Seyfert Is such as MCG-6-30-15 (factor of 3; | Vaughan & 
Fabian|[2004l ) or NGC 4051 (factor of 8 over longer time scales; |Uttley et al.||2004D . This hmited 
flux range prevents strong constraints on the derived value of a, as is clear from the results shown 
in Table [B 

Under the assumption that the spectral variability in NGC 3783 can be described by the sum 
of a constant and variable component where the variable component has a constant shape, i.e. 
a = 1 in equation 1, we can use flux-flux plots for a variety of energy bands with respect to a 
single reference band (Fig. |4]) to obtain the spectral shape of the constant component. Since the 
best fits are indeed roughly consistent with a = 1, this assumption is probably valid. [Taylor et aL 



(2003) and Vaughan & Fabian (2004) describes the methodology used here in detail. 



On inspection, the reference band, F(E1), is chosen to be the 1.0-1.5 keV energy range as this 
shows the least contribution from the constant component (see Fig. [5] and Fig. 14). By measuring 
the offset of the flux-flux plot as a function of the discrete energy band c{E), we can estimate 
the fractional contribution of the constant component to the spectrum as a function of energy. 
Fig. [5] shows the offset c as a function of energy. The offset was normalised relative to the average 
flux in each specific energy band. This normalisation is applied so as to show the shape of the 
constant component independent of the instrument response. Note that by virtue of choosing the 
1.0-1.5 keV as the reference band, we are assuming that this band contributes little to the constant 



-10- 



100 




J I I I I ill III 

1 10 
Energy (keV) 



Fig. 5. — Shape of the constant component as obtained from the hnear flux-flux relations shown 
in Fig. |4} The constant vertical offsets were divided by the flux in each energy bin so as to show 
the shape of the constant component independent of the instrumental response. The vertical line 
shows the expected position of the iron Ka emission hne at 6.4 keV (rest frame). 



component, i.e c(1.0 — 1.5 keV) ^ 0. If, on the other hand, there is indeed a non-negligible 
contribution from the constant component to the overall spectrum at this energy range, then we 
should expect the scale to increase. However, the overall shape of the constant component will 
remain the same. 

From Fig. |5]it can be seen that the constant component is significantly harder than the total 
spectrum in the 2.0-10.0 keV band. A similar analysis for MCG-6-30-15 ( [Vaughan & Fabian 



2004[ see their Fig. 11) shows that the constant component in MCG-6-30-15 has an increasingly 



larger amplitude as the energy softens below 1 keV. This is not the case for NGC 3783 where 
the constant component seems to flatten in this energy range. Furthermore we note the presence 
of a feature at around 6 keV (indicated by the sohd vertical line in Fig. |5]) in the spectrum of the 
constant component as well as a further flattening of the spectrum above lOkeV. These features 



are qualitatively very similar to that found for MCG-6-30-15 by Taylor et al. (2003; see their 



Fig. 5) where the authors identified the constant component as reflection from the inner accretion 



disk, similar to what was proposed by Fabian & Vaughan (2003). 



It should be stressed that the "constant component" identified via this model independent 
method and shown in Fig.[5| needs not be a single component nor does it need to be truly constant. 
It is plausible that this component be in fact the superposition of two distinct components whose 
variability are significantly weaker than that of the "variable" component. 
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Fig. 6. — 0.3-12 keV light curve (1 ks bin) showing the various flux bins used in the following 
section. 



As shown in Paper 1, the time-averaged spectrum of NGC 3783 displays a fairly strong nar- 
row 6.4 keV emission line; this is a strong indication for X-ray reflection from cold and distant 
material. While the exact location of this material is unclear, the limits on the width of the line 
(FWHM < 5500 km s~^; |Shu et al.|2010 ) under the assumption of dynamical equilibrium, place 
this cold material at a distance of r > lO^r^ from the black hole. For a mass of M = 3 x 10^ 
(Vestergaard & Peterson (2006) quotes values ranging from M = (2.8 — 3.0) x 10^ Mq based 
on optical single-epoch spectroscopy), this places the cold reflector at r > 4.4 x 10^^ cm which 
corresponds to 17 light days or more. Thus, during the span of this observation, it is expected that 
the distant/cold reflection component will be approximately constant. Some, but not all, of the 
constant component revealed in Fig. [5] can be identified with this distant reflection. 



3.3. Flux resolved spectroscopy 

A further way to explore the X-ray variability is by examining the evolution of the energy 
spectra with flux. To do this we divided the data into five flux intervals (F1-F5 shown in Fig. [6]) 
each having approximately the same number of photons 2.2 x 10^ and 1.4 x 10^ for the BI 
and FI respectively). Spectra and responses were produced for each flux interval. To illustrate the 
spectral change as a function of flux we show in Fig. [7] both the high (F5; black) and low (Fl; 
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Fig. 7. — Ratio of high (F5; black) and low (Fl; red) flux BI spectrum to an absorbed (F = 1.8) 
powerlaw. The shape of the soft energy range is remarkably similar in both cases suggesting that 
the warm absorber is not changing. As is usually seen in other sources (e.g MCG-6-30-15) the 
overall spectrum becomes harder with a more prominent iron feature at low fluxes. This is as 



expected from the two-component model described in § 3.2 



red) spectra as ratio to an absorbed (A^h= 9.91 x 10^^ cm"^) powerlaw having F = 1.8. It is clear 
that at low flux the spectrum is harder above 2 keV which is as expected from the two-component 
model described in the previous section. Assuming the spectrum is composed of a soft-variable 
component together with a hard-constant component (similar to that shown in Fig. |5]) then as the 
total flux decreases, the relative contribution of the harder component increases. This behaviour is 
also seen in other Seyfert galaxies such as MCG-6-30-15 ( [Vaughan & Fabian|2004[ and references 
therein). 



3.3.1. Difference spectra 

In this section we use the spectra obtained from the five flux levels described above to produce 
flux-resolved "difference- spectra". By subtracting the lowest flux spectrum (Fl) from the remain- 
ing spectra we effectively remove the contribution of any constant (emission) component. This 
technique was described in detail in [Fabian & Vaughan] ( |2003 ) and results in a spectrum consisting 
of the variable component modified by any intervening (Galactic and intrinsic) absorption. 

In subtracting the XIS spectra we assume that the background for each period is constant and 
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Fig. 8. — Left: Data/model ratio of the flux-resolved difference-spectra to an absorbed powerlaw. 
The difference- spectra shown are (from top to bottom) F5, F4, F3, and F2 minus Fl. The spectra 
were fit above 3 keV with a powerlaw and extended to the full energy range after the inclusion of 
neutral Galactic absorption. BI, FI and PIN data are shown in black, red and green respectively. 
Right: Now with the inclusion of a three-zone warm absorber. 



therefore use the difference- spectrum without background subtraction. To test this assumption we 



investigated the background in the XIS-FI data over the five periods described in § |3.4[ In each 
period, the 0.5-10 background count rate are found to be approximately 3 x 10"^ count s"^, which, 
in all cases, corresponds to less that 1.5 per cent of the total source rate. Similarly, the 8-10 keV 
background count rate, although slightly more variable (ranging from 2.3 x 10~^ count to 
2.7 X 10"^ count s"^), always remains at less than 5 per cent of the total source counts. However, 
due to the strong time dependency of the PIN background, we produced background- subtracted 
spectra for each period and then created the difference- spectra using MATHPHA, propagating the 
errors in each arithmetic operation. Figure [8] (Left) shows the data/model ratio for the four in- 
dependent flux-resolved difference- spectra fitted with an absorbed powerlaw above 3 keV. Such 
an absorbed powerlaw results in a satisfactory fit with x^/^ = 1716.5/1669 (1.03) above 3 keV, 
with r = 1.86 ± 0.05 at the 90 per cent level of confidence. Allowing F to vary between the 
various difference- spectra did not result in a statistically significant improvement to the overall fit, 
with a decrease in of 7.8 for three degrees of freedom (F-test value of 2.5). Table |2] details the 
values found for F for the four difference- spectra as well as the global value shown above. When 
we extend the data to cover the full range, we see clear flux deficit below 2 keV due to the warm 
absorber, which is not removed in the difference-spectra. 
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Table 2: Summary of fits to the various flux-resolved difference- spectra. 





Flux-Resolved Difference Spectra 


Flux Bin" 


>3keV^ FuUrange*^ 




Gamma Untied 


Gamma Tied Gamma Tied Gamma Tied 


F5-F1 


1.8 ±0.1 


1.86 ±0.05 2.00 ±0.06 — 


F4-F1 


1.95 ±0.1 




F3-F1 


1.92 ±0.13 




F2-F1 


2.0 ±0.2 








— — 1.94 ± 0.05 




1708.7/1666 (1.03) 


1716.5/1669 (1.03) 2995.3/2853 (1.05) 753.6/709(1.06) 



Notes.- ^ The flux bins refers to the various flux slices shown in Fig. [6j^ The data were fit with a simple 
absorbed powerlaw in the 3-40 keV energy range. ^The full 0.5-40 keV energy range was fit after the 
inclusion of a three-zone warm absorber. ^ T and B refers to flux bins F5 + F4 and Fl + F2 respectively. All 
errors refers to the 90 per cent confidence range. 



In order to model the warm absorber we use the same photoionisation model as in Paper 1 
(generated using XSTAR v2.2.0). The free parameters of this absorption model are the column 
density and ionisation parameter of the absorbing medium. Following from the results of the time- 
averaged spectrum presented in Paper 1 we have modelled the warm absorber assuming three dis- 
tinct zones. To reduce computational time we have frozen the ionization parameter of the various 
zones to the values found in Paper 1. We have also investigated varying the ionisation parameters 
whilst fixing the column densities to the values reported in Paper 1, however it was found to give 
similar results and so we chose to carry on with the former approach. Because the Suzaku/XlS 
detectors do not have the spectral resolution capable of constraining the outflow velocities of the 
various warm absorbers zones, we follow the practice employed in Paper 1 and hold the redshifts 
of these components fixed at the cosmological value for NGC 3783. As such, in each period the 
free parameters of the model were the normalization of the powerlaw together with the column 
densities for the three absorbing zones. Figure [8] (Right) shows the data/model ratio for the var- 
ious difference- spectra after the inclusion of the three-zone warm absorber. The best fit resulted 
in x^/^ = 2995.3/2853 (1.05). Figure [9] shows the fit to the difference- spectra resulting from the 
subtraction of the bottom two flux bins from the top two flux bins (see Fig. [6]). A simple powerlaw 
(F = 1.94 ± 0.05) modified by a three-zone warm absorber again results in an excellent fit to the 
difference- spectra (x^/i^ = 753.6/709 (1.06)). This fit to the full energy spectrum is also sum- 
marised in Table [2} The analysis presented in this section re-enforces the result of the flux-flux 
analysis (§ |3.2| ) that the shape of the variable component (the powerlaw in this case) does not vary 
significantly with flux. 
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Fig. 9. — Similar model to Fig. [8] (right) fitted to the difference- spectra resulting from the subtrac- 
tion of the bottom two flux bins from the top two flux bins (see Fig. |6]). 



3.4. Time resolved spectroscopy 

In order to further visualise the variability we divided the full observation into seventeen 20 ks 
segments (real time) and extracted spectra and responses for each time segment. Each spectrum 
was then compared to the time-averaged spectrum by taking the ratio between them. The result was 
the ratio - at all energies - between the time-averaged spectrum and each 20 ks segment. Figure [TO 
shows this ratio as a function of time and energy (Top), together with the total count rate for each 
segment (Bottom). Apparent changes to the overall spectrum are seen at approximately 60, 180 and 
260 ks, with some indications of a further discontinuity at 300 ks. We note here that from Fig. [10 
it appears possible that there are narrow regions showing spectral variability at timescales shorter 
than 20 ks (e.g. halfway through the fourth period the colour code between 5-7 keV changed 
from orange to green implying a significance greater than 3a), however, spectra extracted at these 
timescales have very low signal-to-noise which prevents us from performing meaningful spectral 
decomposition. We stress that this image is employed only to show the general trend in the energy- 
time parameter space. Based on the broad trends apparent in figure [10] and the light curve shown 
in Fig. [2] we break the full observation into six distinct spectral periods (PI to P6) in order to carry 



out a similar analysis to that of § |3.3| and produce difference- spectra for the five time periods with 
respect to the period with the lowest count-rate (P2 in Fig. [2]). 
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Fig. 10. — Smoothed (left) and un-smoothed (right) time-resolved spectra for 20 ks segments. 
Top: 2-dimension image showing the ratio of the various spectra to the time-averaged spectrum 
as a function of energy and time. The typical standard deviation is approximately 0.03. Bottom: 
Corresponding XIS count-rate. The vertical dotted lines show the position of possible changes in 
the spectral shape as compared to the average. Based on this image and Fig. [2] we divide the total 
observation into six periods. 



3.4.1. Difference spectra 



Figure [TT] (Left) shows the data/model ratio for the five time-resolved difference- spectra fitted 
with an absorbed powerlaw above 3 keV. For the top two difference- spectra (PI and P3 minus P2), 
a powerlaw with an index fixed 1.95, consistent with the value found in the previous section, was 
satisfactory, with x^/^ = 686.0/726 (0.94). Allowing F to vary resulted in slight improvement 
to the fit (Ax^ = 15.6) and a lower value for the powerlaw index, with F = 1.80 ± 0.06. When 
we add the latter three difference-spectra (P4, P5 and P6 minus P2), we find that, whereas the 
3-10 keV range is well fit with a powerlaw having an index of ^ 1.95, the flux above lOkeV is 
severely under-predicted. This is highlighted in the bottom two panels of Fig. [TT] (Left) where 
r = 1.95. Allowing F to vary in the second half of the observation resulted in a statistically 
acceptable fit (x^/i^ = 1793.7/1747(1.03)) however it required uncharacteristically low values for 
F of 1.38 — 1.44. Indeed, if we replace the single powerlaw with a broken-power model having a 
global gamma below 10 keV (F<iokev) and allow F>iokeV to vary between the periods, the various 
difference spectra, including the hard excesses, can be phenomenologically modelled (x^/iy = 
1687.5/1744) with F<iokev = 1.88^^? and F>iokev = 1-8 ± 0.2, 1.9 ± 0.3, 1.071^;}^, 0.731^;}^ 
and O.Srtofs for PI, P3, P4, P5 and P6 minus P2 respectively. It is clear that, whereas the first half 
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Fig. 11. — Left: Data/model ratio for the difference-spectra between (from top to bottom) PI, P3, 
P4, P5 and P6 minus P2. The difference spectra shown on the left are fit between 3 and 10 keV with 
a powerlaw and extended to the full energy range. Right: After the inclusion of a 3-zone warm 
absorber fitted between 0.5-10 keV. In both cases, the P4, P5 and P6 difference- spectra show a 
clear hard excess which is not an artefact of modelling the ionized absorber. BI, FI and PIN data 
are shown in black, red and green respectively. 



of the observation does not require a break in the powerlaw, the second half requires the presence 
of an extra hard component. 



Figure 1 1 (Right) shows the data/model ratio for the various difference- spectra after the in- 



clusion of a three-zone warm absorber similar to that of the previous section (§ |3.3| ). The fit was 
obtained by modelling simultaneously all spectra in the XIS energy range (0.5-10 keV) together 
with the PIN data for period 1 and 3 only. The PIN data for periods 4, 5 and 6 emphasise the level 
of excess above the unbroken powerlaw continuum. The best fit model consisting of a (F = 1.95) 
powerlaw fit to the full energy range but excluding the PIN data for the second half of the observa- 
tion, resulted in xV^ = 3198.7/3189 (1.00). Including the PIN data for periods 4, 5 and 6 results 
in a worst fit (x^/^ = 3518.5/3228 (1.09)) with the excess above lOkeV being the prominent 
source of residuals. Allowing the powerlaw indices to vary between the various difference- spectra 
did not remove the hard excess in the latter half of the observation. We show in Table [3] the values 
obtained for F as a function of time using this simple absorbed powerlaw model. We emphasize 
that since there is still a distinct hard excess in the latter half of the observation, these values are 
only illustrative. 
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Table 3: Summary of fits to the various time-resolved difference- spectra. 





Time-Resolved Difference-Spectra" 


Period'' 


Gamma Untied 


Gamma Tied 


PI -P2 


1.95 ±0.02 


1.95 ±0.05 


P3 -P2 
P4-P2 
P5 -P2 


1 05+0.02 

^•'-'-'--0.02 
-•■•0^-0.03 




P6-P2 


1 07+0. 06 
-'■•O'-0.02 






3496.7/3224 


3518.5/3228 



Notes. -^All cases refer to fits performed in the full 0.5-40 keVenergy range. ^ The difference- spectra are 
made from the periods highlighted in Figs. [2]and[T0j This model successfully accounts for the PIN data 
during the first half of the observation but severely under-predicts the high energy flux in the difference- 
spectra in the latter half. 

The presence of such an excess at hard energies (> 10.0 keV) in the difference-spectrum 
implies that there must be a further varying component w^hose contribution is greatest at these 
energies. Furthermore, this variation does not occur monotonically v^ith flux, since such hard 
excesses w^ere not present in the various difference- spectra presented in § |3.3[ The fact that in all 
difference- spectra the energy range below^ lOkeV w^as w^ell modelled v^ith a simple pow^erlaw^ 
further suggests that this extra component contributes little to the variability below^ this energy. 
A possible origin for the spectral variation at high energies as a function of time could be due to 
change in the ionisation state of the innermost accretion disk and/or change in the contribution to 
the spectrum from the reflection component. The latter can be parametrised by the ratio IZ betv^een 
the reflection and illuminating povv^erlaw^ flux. This possibility Will be investigated further in the 
foUov^ing section. 

3.5. Fitting the individual time segments 

We have shovv^n in the previous sections that the difference-spectra — a rough indication of the 
variable components in the total observed spectrum — can be successfully modelled using a simple 
pov^erlav^ above 3 keV during all flux levels, hovv^ever, w^hen the spectrum is broken dosNXv in time, 
this simple povv^erlavv^ does not describe the latter half of the observation. The difference- spectrum 
for the second half of the observation requires a further variable component which imprints itself 
most notably at hard energies. In all cases, w^hen the pov^erlaw^ continuum is extended to low^er 
energies it reveals the strong presence of a w^arm absorber v^hose ionization and column density 
does not vary w^ith time or flux. A self-consistent origin for both the constant hard-component as 
vv^ell as the varying component above 10 keV (other than the varying pov^erlav^) lies in the intrinsic 
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shape of a reflection continuum arising from the innermost regions around a black hole for a variety 
of ionisations. In the following we adopt such a reflection component as the origin of both the hard- 
constant component and the further variability over the powerlaw component above lOkeV. We 
note that a reflection continuum as such is indeed usually invoked to explain the hard component 
in various Seyfert Is AGN (see e.g. [Zoghbi et aL|2010t [Walton et aL|20T0l ). 



3.5.1. Three-component model 

In order to model the individual time segments as shown in Fig. |2] we construct a model 
based on the primary source of the intrinsic powerlaw continuum being located in a region of 
strong gravity, and gravitational light bending focusing some of the emission from this component 
down onto the disk. The reflection component was calculated with the REFLIONX code by jRoss| 
& Fabian] P005D , convolved with the KDBLUR kernel ( Laor||199T] ) to account for the relativistic 



effects present close to the black hole. In Paper 1 we used the highly sophisticated variable-spin 
relativistic smearing model RELCONV ( [Dauser et al.||2010] ) as the goal was to constrain the spin 
parameter. However, for this study of broadband X-ray variability, such a complex convolution 
model results in very large computational time with no extra benefit beyond that which KDBLUR 
can provide. 

We assume a radial dependence for the emissivity index and ionisation parameter of the disk, 
such that within rbr the disk is highly ionised with a variable ionisation parameter resulting in re- 
flection which is modelled with an inner, hot, relativistically-convolved REFLIONX, and beyond 
this radius it has a constant value of = 1 erg cm s"^ up to an outer radius fixed at 400rg (mod- 
elled with an outer, cold, blurred REFLIONX). After several trials, we fixed the break radius at 
30rg. Changing this value to anywhere greater than approximately 6rg does not change the results 
presented below. To account for the narrow Fe XXVI emission fine present in the spectra we add a 
Gaussian centered at 6.97 keV with a = 10 eV. The normahsation of this line is kept tied between 
each period and we find its equivalent width to vary between approximately 27 eV to 36 eV, in close 
agreement to the value found in the time average spectrum (VI^Fexxvi = 22 ± 5 eV; Paper 1). Re- 
flection from distant, low-velocity, cold matter is modelled with PEXRAV ( [Magdziarz & Zdziarski 



1995| ) together with a further narrow Fe-Ka Gaussian at 6.4 keV. When employing this model, the 
reflection fraction and folding energy of the cold reflector were frozen at T^coid = 0.5 and 200 keV 
respectively, similar to the values found in Paper 1 . We note here that this model differs slightly 



from that employed in Paper 1 in that the previous paper used a superior model (PEXMON; Nandra 



et al.||2007| ) which self-consistently accounts for the Compton backscattered reflection continuum 



as well as the Ka and K/3 emission lines of iron, the Compton shoulder of the iron Ka line, and the 
Ka line of nickel. These fine details are extremely important when modelling the narrow energy 
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range around the broad iron line and thus determining spin, however for the broadband variabihty 
we are interested in this work, the commonly used combination of PEXRAV+GAUSSIAN suffices. 
In all fits, the equivalent width of the iron Ka line lies in the 80-100 eV range (compared with 
= 98 ± 5 eV and 88 ± 6 eV for Suzaku and Chandra data respectively; Paper 1). 

Paper 1 investigated the effect of: (1) fixing the iron abundance of the distant reflector (PEXRAV 
in our case) to that of the inner accretion disk and; (2) allowing these abundances to vary indepen- 
dently. The authors find that a high iron abundance (Fe/solar ?^ 4) is statistically preferred, however 
little difference is found between the two scenarios. Here we fix the abundance of the disk reflec- 
tion (both inner and outer REFLIONX) at 4 times solar, similar to the value used in Paper 1 but 
keep the distant, cold reflector fixed at solar. The Galactic neutral hydrogen column density and 
inner disk inclination were fixed at 9.91 x 10^° cm"^ and 23 degrees as found in Paper 1. The 
powerlaw indices of both REFLIONX components are tied to that of the main continuum. The 
emissivity profile beyond rbr is assumed to follow e{r) oc r"^. We account for the effects of the 
three-zone warm absorber by the use of the XSTAR grids described in § |3.3| with the values for 
the ionisations parameters frozen to that found in Paper 1. Similarly to Paper 1, we also allow 
a fraction / of the continuum to scatter/leak around the warm absorber. This value has consis- 
tently been found to be around 15 per cent. Emission from a black-body with temperature of 
55 eV is also included in the fit. The total model thus consists of seven global parameters - the 
inner disk radius together with the inner disk emissivity index, the powerlaw index F, black-body 
normalization, normalizations of the two narrow emission lines (Fe-Ka and Fe XXVI) and the 
fraction of scattered continuum /. In addition to these seven global parameters, we also have a 
further seven parameters individual to each period - column densities of the three XSTAR grids, 
inner reflection ionization parameter, powerlaw normalization and the normalization of both inner 
and outer reflection components. This model resulted in a good fit to all the individual periods 
{X^/y = 4494.7/3372). 

Figure [12] shows the data/model ratio for all six periods and Table |4] summarises the values 
found for the various parameters of interest. In order to obtain the errors (90 per cent confidence) 
for the various parameters of interest in this study — mainly the inner-disk ionisation, fluxes of the 
various components and the global powerlaw index — we froze the column densities of the warm 
absorber to the best-fit values obtained as we have already shown that the absorption does not 
vary on the time scales probed in this observation. Freezing these parameters was done in order to 
reduce computational time and could result in a slight underestimate of the true error. Nevertheless, 
the global powerlaw index is found to be F = 1.82 ± 0.01, very similar to that found for the time 
averaged spectrum presented in Paper I. We investigated the possibility of a varying powerlaw 
index and show in Table |5] that gamma is perfectly consistent with the global value found above 
throughout the six periods. This supports the results presented in the previous section in which the 
difference spectra below 10 keV is shown to be well modelled with a simple powerlaw. 
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Fig. 12. — Data/Model ratio for the six periods fitted simultaneously with a model consisting 
of a powerlaw, a cold = 1 erg cm s"^) reflection component originating in the outer parts of 
the accretion disk and a "hot" reflection component with variable ionisation between ^100 - 
500 erg cm s"^ shown in the figure) originating from the innermost regions of the accretion disk. 
The continuum is absorbed by the three-zone warm absorber described previously, and the fit 
includes two narrow emission lines (Fe-Ka and FeXXVI) as well as reflection from cold, low 



velocity matter as described fully in § |3.5.1[ See Fig. [14] for a visual depiction of the evolution of 
these three components during the different periods. 



The flux evolution of the three variable components is shown in Fig. [13] (left). It is clear that 
as the [10"^, 10^] keV flux in the powerlaw continuum (blue) begins to decrease, the inner (red) 
reflection also decreases. However, in the latter half of the observation the flux in the powerlaw 
decreases, whereas the inner-reflection flux shows a gradual increase. The reflection arising from 
the outer part of the disk appears to be nearly constant. In fact, an equally good fit (x^/^ = 
4514.6/3395(1.33); i.e. Ax^ = 16.9 for 5 d.o.f) can be achieved if we force it to be constant 
throughout all six periods . The evolution of these components manifests itself in the observed 
evolution of the reflection ratio Tl, where we see in the right panels (top and middle) that the total 
contribution to the spectrum from the reflection component reaches a maximum in the second half 
of the observation and has a clear minimum during Period 2. We also show in the the bottom- 
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Table 4: Summary of fit to the various time-resolved spectra. 
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Notes.- Fit to the individual periods highlighted in Figs.[2]and[T0j The six periods were fit simultaneously 



with a number of global, as well as period-dependent, parameters. The model and the various parameters are 
described in section 3.5.1. ^ The hydrogen column density, disk temperature, inclination and iron abundance 
are frozen (f) at the values quoted in Paper 1. ^ The various fluxes are found between 0.001 and 1000 keV 
and are in units of x 10~^^ erg cm~^ s~^. 

Table 5: Summary of fits having the power law both tied and untied between the six periods. 



Gamma 


Period 1 


Period 2 


Period 3 Period 4 


Period 5 


Period 6 


X^/d.o.f. 


Tied 
Untied 


1.82 ±0.01 


-1 09+0.02 

i.oz_o 03 


1.82 ±0.01 

1 oq+0.02 1 oq+0.02 

i.OO_o 01 -'-•^^-0.01 


-1 09+0.02 
-•-•^^-0.01 


1.83 ±0.02 


4497.7/3390 
4487.4/3385 



right panel the evolution of the inner ionisation parameter = Lint/nR'^). Assuming the intrinsic 
luminosity (Lint) of the powerlaw continuum (as seen by the accretion disk) as well as the hydrogen 
number density of the disk (n) remains constant then the maximum change in ionisation observed 
{ihighl^iow ^ 4) implies an increase in the height of the corona (above the disk) by a factor of 
^ 2 (assuming Newtonian geometry). Of course it is likely that there could be some contribution 
from both intrinsic variations in the continuum luminosity, and possibly changes in the hydrogen 
number density of the disk as well. In Figure [14] we visually depict the evolution of the three 
components during the course of the observation. The greater contribution from the inner reflection 
component (red) to the hard energy band during the last three periods in comparison to the first 
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Fig. 13. — Left: Flux evolution of the three variable components with the powerlaw, inner and 
outer reflection shown in blue, red and green respectively. Right: Ratio of the inner reflection (top) 
and total reflection (middle) to the powerlaw flux. The bottom panel shows the evolution of the 
ionisation parameter ^, 
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Fig. 14. — Best- fit model for the six periods, as described in ^ 3.5.1 The total model is shown in 
black, with the powerlaw, inner and outer reflection components shown in blue, red and green re- 
spectively. The global blackbody, absorption and narrow emission lines were removed for display. 
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three is clearly seen (the near constancy of the reflection component arising from the outer parts 
of the accretion disk - shown in green - can be used as a reference point to highlight the level 
of change in the inner reflection flux). The extra flux at hard energies coming from the increase 
of the inner reflection component in comparison to the powerlaw continuum directly explains the 



presence of a hard excess in the difference- spectra presented in §3.4| for the second half of the 
observation in comparison to Period 2. 

4. Discussion and Conclusions 

The main results concerning the spectral variability of NGC 3783 as determined from the long 
Suzaku observation can be summarized as follows: 

1 . The spectral variability of NGC 3783 can be described by the superposition of a soft- variable 
component together with a quasi-constant hard component. The shape of this quasi-constant 
component (Fig. [5]) is hard in the 3-10 keV range as found for several other AGN, and flat 
both below 1 keV and above 10 keV. These features are qualitatively similar to that expected 
from reflection arising from the inner parts of an accretion disk. 

2. A simple powerlaw provides a good fit to all the difference-spectra between 3-10 keV sug- 
gesting that in this band the variability is dominated by variation in the normalization of the 
powerlaw-like Comptonized continuum. 

3. Extending the powerlaw to lower energies reveals the strong presence of a warm absorber 
(Figs. |l|8| and[TT]). The warm absorber does not vary with time or flux (see e.g. FigjT]), which 



is consistent with previous studies of this source. 

4. The presence of complex spectral evolution is clearly shown in Fig. [10} The evolution of the 
softness ratio allows for the clear distinction of at least six "spectral periods" (Fig. [2]). 

5. Time-resolved difference- spectra show the presence of an excess flux to a powerlaw contin- 



uum above lOkeV in the latter half of the observation (Fig. 11 ). This excess is not seen in 
the flux-resolved difference- spectra (Fig. [8]) suggesting that the further variable component 
above 10 keV varies with time but not monotonically with flux. 

Changes in the reflection fraction and possibly ionisation state of the inner accretion disk 
could result in a reflection component having a quasi-constant level below lOkeV but dif- 
ferent contributions at higher energies . A model consisting of such a reflection component 
together with a powerlaw and the associated Compton backscattered emission from cold ma- 
terial, provides an excellent description of the continuum for the six distinct spectral periods 
seen in NGC 3783 (Fig. [12]). 
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Our time-domain spectral analysis presented in Section |3.5| presents a complex picture. To 
organize this discussion, we choose to focus on three primary observables, the observed power- 
law flux (PL), the ratio of the total disk reflection to the powerlaw flux TZ = {Refl + Ref2) / PL, 
and the ionization state of the inner accretion disk (Q. Thus, we have a 3-dimensional parameter 
space (PL, TZ, ^) that, broadly, describes the X-ray flux and spectral state of the central engine. 

What are the theoretical expectations for the evolution of the source in this parameter space? 
We suppose for now that the observed continuum X-rays are emitted from a single structure ("the 
corona"), and that the observed reflection spectrum results from irradiation of the disk by that same 
structure. Within this picture, we can imagine at least three distinct causes for flux- and spectral- 
variability. Firstly, suppose that the X-ray emitting corona increases its luminosity while all other 
physical properties of the system (coronal geometry, coronal beaming, and disk density) remain 
constant. Clearly both the observed flux and the ionization parameter of the inner disk would in- 
crease in proportion, i.e., we would expect ^ oc PL. Since the geometry is unchanged, however, 
we would expect TZ to remain constant. The second type of variability corresponds to a change 
in the location, geometry, size or beaming pattern of the corona while its luminosity and the disk 
properties remain unchanged or changes in the geometry of the disk in a manner that some discrete 
Compton-thick clump of gas could move into the irradiation zone thus changing the (anisotropic) 
scattered/reflected emission directed into the line of sight. In the case that these changes lead to 
increased irradiation of the disk at the expense of flux reaching the observer (through, for exam- 
ple, increased gravitational light bending if the source moves closer to the black hole), PL would 
decrease while both ^ and TZ would increase in rough proportion. The final type of change that 
we consider are density changes in the surface layers of the accretion disk (with all other aspects 
of the system remaining unchanged). Such behaviour may result from viscous time-scale evolu- 
tion^ following some dramatic accretion or ejection event. In this type of change, ^ would respond 
inversely to density changes, while TZ and PL would remain unchanged. 



Figure [B] shows the evolution of NGC 3783 projected onto the (PL, ^), (PL, 7^), and (7^, ^) 
planes. Clearly, none of the three simple scenarios presented above can describe the complex 
behaviour that we see — the variability must be due to a combination of physical changes. The 
fact that TZ changes demonstrates that there must be structural and/or geometrical changes in the 
X-ray emitting corona and/or accretion disk. The most dramatic of these changes are in the final 
three data segments, although there is some weaker (and possibly flux-correlated) variability of TZ 



even between the first three data segments (Fig. 15 , right panel). 



VISCOUS 



^For NGC 3783 the dynamical time-scale ^) at 6rg is approximately 350 s. Assuming 
{H / R)~'^Vt~^ and a = (H/R) = 0.1 we have Tyiscous ^ 350 ks, similar to the duration of the observation 
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Fig. 15. — Variability of NGC 3783 in the parameter space defined by the observed power-law 
flux PL (in units of xlO~^°erg cm"^ s"^, ionization parameter of the inner disk ^ (in units of 
erg cm s"^), and the ratio of the total reflected flux to the observed power-law TZ = {Refl + 
Ref2)/PL. Lines connect data points from consecutive time segments; the first time-segment has 
a filled marker. In the left panel, the dashed line shows a representative line with ^ oc PL expected 
from pure luminosity evolution of the X-ray corona; clearly this is a poor match to the data. The 
dotted fine shows the best-fitting relationship ^ oc {PL — PLq) where PLq is some constant flux 
component that has little or no associated reflection. 

It is interesting to note that ^ and PL do indeed seem correlated. However, we can only 
recover the linear proportionality between ^ and PL (expected for "pure" coronal luminosity 
changes) if we assume that there is a constant flux component (corresponding to PL = 2.63 x 
2Q-10 (3321-2 g-i^ ^j^^^ j^^^ ^^^y weak associated reflection spectrum, e.g., highly beamed 
emission from a steady jet. In this picture, all of the observed flux variability is due to the corona 
which also is the principal source irradiating the disk. In addition to reproducing the observed 
^-PL correlations, it also explains the weak correlation of TZ with PL during the first three time 
segments — as PL increases, both its contribution and the contribution of its associated reflection 
spectrum increase relative to the constant component. The dramatic increase in TZ at the end of the 
dataset still requires structural changes, however. 

Like any data, the diagnostic power of these spectra are necessarily limited. However, a rela- 
tively simple and physically-motivated analysis of the spectrum and source variability shows that 
reflection from a disk with corresponding changes in ionization provides an excellent description 
of the data. This discussion illustrates the power and utility of long X-ray datasets for examining 
AGN spectral variability. Such studies have the potential to yield new insights into the complex 
interplay between various physical processes in the cores of active galaxies. 
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